Critical state in thin anisotropic superconductors of arbitrary shape 
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A thin flat superconductor of arbitrary shape and with arbitrary in-plane and out-of-plane 
anisotropy of flux-line pinning is considered, in an external magnetic field normal to its plane. It 
is shown that the general three-dimensional critical state problem for this superconductor reduces 
to the two-dimensional problem of an infinitely thin sample of the same shape but with a modified 
induction dependence of the critical sheet current. The methods of solving the latter problem are 
well-known. This finding thus enables one to study the critical states in realistic samples of high-T c 
superconductors with various types of anisotropic flux-line pinning. As examples, we investigate the 
critical states of long strips and rectangular platelets of high-T c superconductors with pinning either 
by the ab planes or by extended defects aligned with the c axis. 



PACS numbers: 74.60.-w, 74.60.Ge, 74.60.Jg 



I. INTRODUCTION 

Numerous publications deal with the critical state 
problem in isotropic thin superconductors placed in a 
perpendicular magnetic field.EI By thin we mean that 
the sample thickness d is considerably less than its lat- 
eral extension. Analytic solutions of the problem were 
obtained in the cases of a circular diskcl, a stripO, and 
elliptic-shaped films .□ For thin superconductors of an ar- 
bitrary shape, numerical methods for solving this prob- 
lem were elaborated in Ref. |],[| and the critical state of 
a rectangular platelet was investigated in detail in Ref. |7]. 
All the above-mentioned results can be derived consider- 
ing the appropriate flat superconductors to be infinitely 
thin. In the framework of this approach a possible 
plane anisotropy of pinning can be taken into accountpEl 
yielding good agreement with appropriate magnetooptic 
experiments!] 

It is important to note that in the highly anisotropic 
high-r c superconductors the flux-line pinning in general 
depends on the angle 9 between the local direction of 
the magnetic induction B and the c axis which in typical 
experiments is normal to the plane of the samples. For 
example, this type of anisotropy occurs when one takes 
into account the intrinsic pinning exerted by the CuO 
planes.u Besides this, twin boundaries, columnar defects, 
and other extended defects give rise to such an out-of- 
plane anisotropyJd In all these cases, the flux-line curva- 
ture that always occurs in thin superconductors placed in 
a perpendicular magnetic field, leads to a dependence of 
the critical current density j c on the coordinate z across 
the thickness of the sample. Therefore, the critical state 
problem becomes three-dimensional (3D), and the feasi- 
bility of treating such superconductors as infinitely thin 
requires special consideration. 

For a thin anisotropic disk with radius R placed in 



a perpendicular magnetic field it was showiS that the 
small ratio d/R enables one to split the 2D critical state 
problem into two one-dimensional problems. The first 
one treats the critical state across the thickness of the 
disk, while the second one treats the disk as infinitely 
thin and gives the distributions of magnetic field and 
sheet current along its radius. Using this method, the 
magnetic moment of a thin anisotropic disk was theoret- 
ically studiedO in the case when the angular dependent 
j c has a peak at 8 = 0. In this special geometry of the cir- 
cular disk the directions of field and current are fixed and 
known in advance, and the critical state problem thus is 
essentially not 3D but only 2D. This simplification does 
not take place in thin samples of arbitrary shape, where 
the current stream lines, magnetic contour lines, and pen- 
etrating flux fronts in general do not coincide.u Moreover, 
according to Ref. [|, in the general case the directions of 
the circulating currents may change across the thickness 
of the sample, and a rotation or twist of the flux line 
arrangement is possible. 



In the present paper we show that even in the general 
3D critical state problem of a thin superconductor with 
arbitrary shape and arbitrary anisotropy of j c , the 3D 
pinning problem can be split into two simpler problems 
which can be solved: The ID problem across the thick- 
ness of the sample is solved analytically, while the in- 
plane problem reduces to that of the infinitely thin super- 
conductor which was treated before. Since the methods 
of solving the latter problem are well elaborated^El our 
finding enables one to study the critical state in flat sam- 
ples of arbitrary shape and with any type of anisotropy. 
As concrete examples, we shall investigate the critical 
states of long strips and rectangular platelets with vari- 
ous types of out-of-plane anisotropy. 
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FIG. 1. The front of magnetic flux penetrating into a thin 
rectangular superconductor plate with pinning when the ap- 
plied perpendicular magnetic field increases from zero. The 
top figure shows the two-dimensional curve 70 that forms the 
equator of the three-dimensional flux front 7 shown in the 
lower plot. Inside this 3D front the magnetic induction B 
and current density j are exactly zero, and inside the curve 
70 the perpendicular component B z is practically zero when 
d <IC (a, b) for a thin plate with dimensions 2a x 26 x d. 



II. SPLITTING THE 3D PROBLEM 



In what follows we assume that the pinning-caused 
characteristic magnetic field H cr = j c d/2 considerably 
exceeds the lower critical field H c \, and the thickness d 
of the sample is much greater than the London pene- 
tration depth A. Under these conditions the magnetic 
induction is practically equal to the magnetic field in the. 
superconductor, and the so-called geometrical barriered 
is negligible. We study the critical state problem in 
the framework of the macroscopic approach, in which all 
considered lengths are larger than the flux-line spacing, 
and treat the superconductor as a uniform anisotropic 
medium. 

Let us place the coordinate system so that its xy plane 
coincides with the middle plane of the sample and the 
z axis is directed along the external magnetic field; the 
case of a rectangular platelet is shown in Fig. 1. Since 
the thickness d of the sample is assumed to be much less 
than its smallest lateral dimension L, the strict equations 
V • j = and V • B — give that, in leading order in the 
parameter d/L, one has j z = 0, and the perpendicular 
field component B z is independent of z inside the sam- 
ple, i.e., B z = B z (x,y). Accounting for the fact that the 
scale of changes with x and y of B z and of the in-plane 
field B t = (B x ,B y ) is L, while the appropriate scale for 
changes of B t across the thickness is d, we may neglect 



the derivatives of B with respect to x and y. The critical 
state equation can thus be written in the form: 

z x -5— = MoJc , 

oz 

where z is the unit vector along z and j c is the criti- 
cal current density. This equation means that the criti- 
cal gradient develops predominantly across the thickness 
of the sample, and we can consider the fields and the 
currents to be independent of x, y in sufficiently small 
parts of the sample with dimensions in the x, y direc- 
tions much greater than d and much less than L. To 
proceed our analysis we write Bf = B t t and j c = j c n, 
where n = (cos<^, sint^) and t = (costp, simp). The 
angles tp(x,y, z) and ip(x,y,z) define the directions of j c 
and B t in xy planes. The above critical state equation 
can then be rewritten as 



dBt ... 

-q^ = MoJ c sin(p - ip) , 

dip 

B t — = -fi j c cos(ip-ip). 



(1) 
(2) 



We emphasize that in the general case the critical cur- 
rent density j c may depend on B z , B t , tp — tp, ip, i.e., 
jc = jc(B z , Bf,tp — ip,ip). The quantities B t and B z can 
be also expressed in terms of the magnitude B = |B| and 
the tilt angle 9 between the local direction of B and the 



B, 



B cost 



B t = BsinO. 



The dependence on ip occurs only if there is an anisotropy 
of pinning in the a;^,plane (e.g., when twin boundaries 
exist in the sampleO). In the isotropic case, when the 
flux-line pinning depends neither on tp nor on the tilt 
angle 9, the dependence of j c on B z and B t occurs only 
through the magnitude B. Note also that, since B z is in- 
dependent of z, it enters Eqs. (1,2) as a parameter, and 
thus we have only two equations for the three functions 
B t , tp, and Tp. One more equation can be obtained if one 
takes into account the prehistory of the critical state, 
telling how the flux lines were pushed into the sample. 

In increasing moderate magnetic field there is a flux 
and current-free core in the superconductor, Fig. 1. The 
surface of the core, 7, is the penetrating flux front, which 
we describe by a function z = z~f(x,y). In the limiting 
case of an infinitely thin superconductor (d — > 0) the 
flux front may be thought of as a flat curve 70, forming 
the outer rim (equator) of 7, since we have B z = in- 
side this contour. In the case of small but finite d <C L, 
the component B z is still practically zero in the region 
inside 70 since there the flux lines are almost parallel 
to the flat surfaces of the sample. Note that these flux 
lines may have different orientations at different values 
of z if the sample is not a circular disk or a strip;cl 
thus, in principle, flux^cutting processes may occur in 
the superconductor .EJ~t£l 

We now consider the third required equation in the re- 
gion Zj(x,y) < \z\ < d/2 with x, y inside the curve 70, 
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i.e. between the lens-shaped empty core and the flat sur- 
faces. In increasing applied field the flux lines are pushed 
continuously from the flat surfaces into this region, such 
that the in-plane orientation of the flux lines is a unique 
function of B t , x, and y: 



ip = V(Bt,x,y). 



(3) 



The function '5 is found from the solution of the critical 
state problem for the infinitely thin sample (see below). 
Equation (3) is just the third required equation. Now we 
are able to solve Eqs. (1-3). Dividing Eq. (2) by Eq. (1) 
we obtain 



dip 

Bt ~dB t = - cot (^-^)• 



(4) 



Inserting here the known field orientation, Eq. (3), we 
find the difference tp — ip. and the angle of the currents tp 
as functions of B t . Finally, the dependence of B t on z 
can be obtained from the implicit relation 



db. 



o tM)jc(0,b t ,Tr/2,ty(b t ,x,y)) ' 



(5) 



which follows from Eq. (1), and from the formula 



j c (0,B t ,ip-ip,ip) 



sin(y> — ip) 



reflecting the fact that only the component of j c normal 
to B is essential (in absence of flux cutting). The function 
Zj(x, y) describing the shape of the flux front is obtained 
from the equality: 



d 
2 



dh 



o Mojc(0,6 t , tt/2, V(b t ,x,y)) ' 



(6) 



where B s (x,y) = B t (x,y,d/2) is the tangential compo- 
nent of the surface field. This B s may be obtained by 
solving the critical state equation for infinitely thin su- 
perconductors. Thus, Eqs. (3-6) give the solution of the 
critical state problem across the thickness of the sample 
in the region inside the curve 70. 

It should be emphasized that the above formula for the 
dependence of j c on tp—ip and Eqs. (3,5,6) hold only if flux 
cutting does not occur in the superconductor when the 
critical state is established in it. Flux cutting will only 
happen if \cot(p — ip)\, calculated from Eq. X4), turns 
out to be greater than some critical value x«iJp this 
case the so-called generalized critical-state modeEj must 
be used to describe the critical state in the region of the 
sample inside 70. However, since at present we cannot 
point out the shape of the superconductor for which this 
flux cutting really occurs, we do not analyze this case 
here. 

Let us now derive the third equation for the region 
where B z ^ 0, i.e., outside the curve 70. In deriving it is 
necessary to take into account that in an anisotropic su- 
perconductor the direction of the flux- line velocity in gen- 
eral does not coincide with the direction of the Lorentz 



force causing this movement, and thus the appropriate 
electric field E is not always parallel to the current den- 
sity (see Appendix A). With this in mind, the third equa- 
tion results from the following two conditions: First, the 
shape of a flux line cannot change when an applied mag- 
netic field H increases slightly so that 70 and all the flux 
lines in this region of the sample are shifted by distances 
which are much less than the sample width L. Second, 
during the shift, each line element moves along the local 
direction of [E x B] . The geometrical analysis of these 
conditions gives 



B 2 Z sin £ — B\ cos ip s'm(ip — £ ) 
B z cos £ + B? sin ip s'm(ip — £) 



= const 



where the angle £ is determined by Eqs. (A2,A3,A5) in 
Appendix A. Since B t (0) = 0, the obtained relation can 
be rewritten as 



tan[£ - £(0)] 



B} su# - £(0)] cos[</> - £(0)] 



Bf cos 2 [ip 



(7) 



with £(0) = £(x, y, 0). This is the third required equa- 
tion in the region where B z ^ 0. When an in-plane 
anisotropy of the flux line pinning is absent, it is pos- 
sible to find the solution of Eqs. (1,2,7) in the interval 
< z < d/2. The initial conditions to these equations 
are 



B t (0) = 0, tp(0) - r/>(Q) = tt/2 . 



(8) 



where the last equality follows from Eq. (2) with B t (0) = 
0. Since the solution of the equations is unique, it is suf- 
ficient to guess it. The solution has the form 



V{z,V,z) = fo(x,y) , 
ip(x, y, z) = (p (x, y) - tt/2 = ip (x, y) , 

dh 

V>ojc{B z ,b t ,ir/2,'ip ) ' 



(9) 
(10) 

(11) 



Here the angle <po(x,y) defines the direction of the sheet 
current in the critical state of the infinitely thin super- 
conductor and is assumed to be known. In obtaining 
the solution we have taken into account that £ = ip at 
ip — ip = tt/2 when the in-plane anisotropy is absent. 
Eqs. (9-11) mean that the flux lines are curved in the 
region where B z ^ 0, but they are not twisted there. In 
other words, the distributions of the magnetic induction 
and the current across the thickness of the sample are 
the same as in a circular disk or a strip if one considers 
small regions in the xy plane with the same values of B z . 
Eq. (11) allows us to find an implicit relation between 
the sheet current J and B z : 



J/2 



dht 



jc(B z ,fJ. h t ,7T/2, ip ) ' 



(12) 



3 



Here we have used that B t (z = d/2) = B s = noJ/2, 
where J = |J| and J is the current density integrated 
over the film thickness. It should be noted that although 
we indicate explicitly the angle tpQ as an argument of j c in 
Eqs. (11,12), the critical current density is independent 
of this angle in cases without in-plane anisotropy. 

We now discuss the case when some in-plane 
anisotropy of flux-line pinning exists in the superconduc- 
tor together with the out-of-plane anisotropy, i.e. we al- 
low for a dependence of j c on ip. This case can occur, e.g. 
in twinned crystals of high- T c superconductors. When an 
in-plane anisotropy exists, it is impossible to find the so- 
lution of Eqs. (1,2,7) in the general form. We may only 
state that in this situation the in-plane angle tp generally 
depends on z, and thus flux lines are not only curved but 
also twisted. Then, according to Eq. (2), the current is 
not normal to the magnetic induction. The magnitude 
of the twist is determined by the relative strengths of the 
in-plane and out-of-plane anisotropies (and, of course, by 
the thickness of the sample) . Interestingly, when j c does 
not depend on B t , the solution described by Eqs. (9-11) 
becomes true again, and the twist disappears. (This sit- 
uation occurs, e.g., in sufficiently large magnetic fields, 
H 3> H cr , or if there is no out-of-plane anisotropy in the 
superconductor and j c is practically independent of the 
magnitude of the magnetic induction). However, if H 
changes, in this case a flux line is shifted in a direction 
which does not lie in the plane containing the line. Ac- 
cording to the results of Appendix A, the direction makes 
an angle a with the above-mentioned plane where 



tana = — 



d[lnj c (B z ,ip - ip ,ipo)] 



dip 



tp = "00+ 



To sum up, it is important to emphasize that at a 
given dependence of j c on B z , B t , <p — tp and tp, the one- 
dimensional equations (1), (2), (7) with initial conditions 
(8) can be solved (at least numerically) in the case of an 
arbitrary anisotropy of flux-line pinning, and a relation 
J(B z ,ipo) generalizing Eq. (12) can be obtained. Here <po 
defines the direction of the sheet current in the xy plane. 

Thus, we have reduced the 3D critical state problem 
to a 2D problem in which the independence of the 
sheet current J(B z ,(po) is determined not only by the 
-^-dependence of the critical current density j c but also 
by its dependence on the tilt angle 9 of the flux lines 
away from the film normal, see Eq. (12). It should be 
noted, however, that to obtain this reduction we have as- 
sumed that d is much less than the characteristic scale of 
changes of B in the xy plane. Numerical solutions of the 
critical state equation for infinitely thin superconductors 
of various shapes show that this assumption indeed is 
true everywhere; the scale of variation is of the order of 
L except for a small region near the flux front 70 where 
the derivatives dB x /dy and 8B y /dx are not small com- 
pared with dYSt/dz. However, the width of this region is 
determined by the thickness d and is thus small. 

Knowing the function J(B z ,ipo) obtained from 



Eq. (12) or from the above-mentioned generalization of 
this equation, it is possible to solve the critical state 
equations., for an infinitely thin superconductor of any 
shape.ErErEI The static equations to be solved are 



V-J = 



(13) 



and the Biot-Savart law connecting B z with the sheet 
current J(x,y) = J_f /2 j(x,y, z) dz, 



^ 1 B z (r)=H 



Is 



[r x j; 

47ri? 3 



■d 2 i 



(14) 



where R = r — r'; r and r' are 2D vectors in the xy plane, 
and the integration is carried out over the specimen area 
S. In addition, one has to take into account that B z = 
inside the 2D flux front 70, while outside 70 the sheet 
current is equal to J{B z ,ipo). At given applied field H, 
a solution of this static problem exists only for a certain 
curve 70, which determines the shape and position of the 
penetrating flux front. In Refs. HJlS m ^ &c ^, the more 
general dynamic equations were considered, which allow 
one to find the solution of the 2D critical state equations 
very effectively for any prehistory H(t). The appropriate 
formulas are presented in Appendix B. 

Solving these equations for the infinitely thin super- 
conductor of arbitrary shape in a perpendicular applied 
field H, one obtains B z (x,y, H), J(x,y,H), ipo(x, y, H), 
B s (x,y, H) — noJ/2, the direction of B s , i.e., the func- 
tion tpo(x, y, H) = (fio(x, y, H)—tt/2, and also the position 
of 7o- Eliminating H from tpo(x, y, H) and B s (x, y, H) in 
the region inside 70 , and replacing B s by B t we arrive at 
the function ^{B t ,x,y) that enters Eq. (3). Thus, the 
solution of the 2D critical state problem for the infinitely 
thin superconductor enables one to obtain all the neces- 
sary information for the calculations of the core shape 
z-y(x,y) and distributions of B t , tp, if across z at any x 
and y. The knowledge of all these functions permits to 
describe the 3D critical state of thin anisotropic super- 
conductors in detail. 



III. OUT-OF-PLANE ANISOTROPY 

We now use the results obtained above to investigate 
the critical states of thin superconductors with out-of- 
plane anisotropy alone. That is, we consider the case 
when the critical current density j c (B z , B t , 7r/2, tp) does 
not depend on the in-plane angle tp. In addition we here 
neglect its dependence on B such that j c depends only 
on the flux-line tilt angle 9, j c — j c {S). This situation 
occurs when the flux-line pinning is isotropic in the xy 
plane, and the scale Bq characterizing the dependence of 
j c on B considerably exceeds the self-fields of the critical 
currents. As shown above, the anisotropy of j c (9) can be 
accounted for by introducing an appropriate B z depen- 
dence of the sheet current J = J(B Z ), implicitly given by 
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FIG. 2. The out-of-plane anisotropy of the critical current 
density j c (8) described by the model (18) (top plot) and the 
equivalent induction dependence of the sheet current J(B Z ) 
of Eq. (17) (lower plot). In both plots the upper half is for 
the parameter p = 1 and the lower half for p — —0.8. The 
second parameter is q — 0.05 . . . 50. 

Eq. (12). In the considered case, a more explicit form of 
this dependence may be obtained in the following way. 
Differentiation of Eq. (12) with respect to B z yields the 
two equations: 



Jc (9)d = J(B z )-B 
tan 6* = 



dJ{B z ) 
Z dB z 

2B Z 



(15) 
(16) 



In fact, Eqs. (15,16) are the parametric form of that func- 
tion j c (0) which results in a given dependence J(B Z ) via 
Eq. (12). 

To fix ideas, we consider the following rather general 
model dependence for J(B Z ): 



J(B z )=j c (0)d 



1+p exp — 



qB z 



(17) 



where B cr — jj,oj c (0)d/2, while j c (0), p and q are the 
parameters of this example (p and q are dimensionless). 
One then easily verifies that the corresponding angular 
dependence of the critical current density takes the form: 

jc(0)=jc(P) [l+p(l + qt)exp(-qt)] , 

tan 9 = r 1 [1 + p exp(-g t)] , (18) 



where t is a curve parameter with range < t < oo 
equivalent to n/2 > 9 > 0. This model dependence j c (9) 
is presented in Fig. 2 together with the corresponding 
J(B Z ), Eq. (17). For appropriate choices of the param- 
eters p and q this model describes the intrinsic pinning 
by the CuO planes in high-T c superconductors (p > 0, 
maximum j c at 9 = tt/2) and pinning by columnar de- 
fects perpendicular to the film (p < 0, maximum j c at 
9 = 0). It may thus be used to simulate these impor- 
tant cases. We remark that the critical state of a circular 
disk with some model for intrinsic pinning was studied 
numerically in Ref. |l7], while the magnetization of the 
disk in the case when j c has a peak at 9 = was ana- 
lyzed in Ref. |ll|. However, in those papers only the fully 
penetrated critical state was considered. 

Below we calculate the partially and fully penetrated 
critical states in strips and rectangular platelets with 
anisotropic pinning using the angular dependence j c {9) 
described by Eqs. (18). The corresponding results for cir- 
cular disks are ver*-similar to the results for the strip. In 
fact the algorithmlil used below to time-integrate the ID 
equation of motion for the sheet currentp/(x, t) in thin 
strips is almost identical to the algorithmta used for thin 
disks. 



IV. THIN STRIP 



We begin with the calculation of the critical state in an 
infinitely thin strip (see Appendix B). In this calculation 
the function J{B Z ) described by Eq. (17) is considered 
as the critical value of the sheet current. The obtained 
results are presented in Figs. 3-7. 

Figure 3 shows the induction profiles B z (x) and sheet 
current profiles J(x) in increasing H for thin strips with 
various anisotropies of j c (9) depending on the param- 
eters p and q. The profiles are presented for several 
values of xq, the position of the 2D flux front 70. The 
reference case of isotropic pinning (p = 0) is also shown 
here. In Fig. 3a one can see that for p > (pinning by 
the CuO layers) the penetrating flux-density profile is 
very steep at x = xq: B z (x) jumps almost abruptly and 
J{x) has a very sharp peak which ideally should reach 
the maximum value J(xq) = (1 + p)j c (0)d, cf. Eq. (17). 
In our computation this ideal maximum height is not 
reached due to our finite creep exponent n = 101 and 
finite number = 300 of equidistant grid points (to 
illustrate the influence of n the case of a small value 
n = 11 is shown in Fig. 3d). As opposed to this, in 
Figs. 3b and 3c the induction profiles for p < (pin- 
ning by linear defects perpendicular to the strip plane) 
are less sharp: at x = xq, B z (x) vanishes less steeply 
(almost linearly for Fig. 3b) and J(x) decreases mono- 
tonically and exhibits an inflection point with infinite 
slope, dJ/dx\ x=Xo = 00, J(xo) = (1 + p)j c (0)d. These 
general features of the flux front of superconductors with 
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FIG. 3. Spatial profiles of the perpendicular field component B z (x) (upper plot) and of the sheet current J(x) (lower plot) 
of a thin strip with anisotropic pinning described by the model (18), see Fig. 2. The various curves correspond to increasing 
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FIG. 4. The x coordinate xq(H) of the penetrating flux 
front, or the penetration width a ~ xq, for thin strips with 
various anisotropies of pinning jc{9), Eqs. (18), calculated 
with large creep exponent n = 101. The line marked by dots 
shows the static analytical result for isotropic strips, Eq. (20). 

anisotropic pipaing are derived analytically in a forth- 
coming paper .E3 It turns out that in the vicinity of xo the 
profiles are B z (x) oc (x— x®) 13 and J(x)—J(xo) <x \x— xo] 13 
where (3 ~ 0.5 — ir~ 1 arctan(pq). 

In the isotropic case (p = 0) our numerical results-prac- 
tically completely coincide with the exact solutionH 



J(x) 2 | 
= — arccot< 

Jed 7T I 



0, 



(4 



x 2 )a 2 - 



(a 2 

Xq = 



x1)x 2 



cosh(7r H/ j c d) 



(19) 
(20) 



where < x < a and 2a is the width of the strip. It 
should be also mentioned that in the case of monotoni- 
cally decreasing J C {B Z ) the critical states of a thin strip 
and a thin disk were studied in Refs. 
In Ref. M the Kim model J C (B Z 



21 2 



respectively. 
J C (Q)B /(B Q + \B Z \) 
was considered, while in Ref. |22| the Kim model and the 
exponential dependence J C (B Z ) = J c (0) exp(— \B Z \/B ) 
were analyzed where J c (0) and Bq are some constant pa- 
rameters. Our results in the case p > are similar to 
the results in-these papers. Equations (15,16) show that 
both modelsEHE3 may be interpreted as angular depen- 
dences of j c {9) with maxima at 9 — tt/2 if Bq is of the 
order of J c (0). 

The position x = x of the flux front 70 as a function 
of the increasing applied field H is shown in Fig. 4 for 
the same anisotropic strips. One realizes that for p > 
the penetration depth a — xq is smaller, and for p < 
larger, than for the isotropic case p — 0. The isotropic 
case is shown twice in Fig. 4: The numerical result for 
xq obtained with a finite creep exponent n = 101 is only 
slightly smaller than the ideal static result, Eq. (20). This 
demonstrates the accuracy of our computations. 

Using the above results for the infinitely thin strip and 
equations of Sec. II, we can now describe the two dimen- 



sional critical state of the anisotropic strip with small 
but finite thickness. 
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FIG. 5. Flux fronts z^(x) in thin strips in increasing field 
H calculated from Eq. (21) and the data J(x) for isotropic 
pining (p — 0, dashed lines) and anisotropic pinning with 
p = 1, q = 4 (solid lines) and p — —0.5, 2 (dash-dotted lines). 
Shown are two applied fields H = 0.1, 0.5 in units of j c (0)d. 
Note the different scales for x and z. 



The general 3D flux front 7 (Fig. 1) is determined by 
Eq. (6) which is rewritten for our model Eqs. (17,18) as 



J(x,y) 



J- 



70 



(21) 



where J 7o = (1 + p)j c (0)d and J(x,y) is the value of the 
sheet current at a given point inside 70. Fig. 5 shows 
z~[(x), the spatial profile or cross section of flux fronts 
7, in thin strips with anisotropic pinning. Note that for 
p > the fronts have a wide flat part near the central 
plane z = 0, while for p < the fronts have a more 
rounded shape near x — xq as compared to the isotropic 
case p = 0. 

Physically, the wide flat flux front for p > results 
from the fact that in this case the flux lines oriented per- 
pendicularly to the film (and to the CuO planes) are less 
pinned and can thus penetrate more easily than the flux 
lines oriented parallel to the film plane. In the opposite 
case p < 0, flux lines in the film plane move more easily 
and cause a more pointed wedge-like front. 

The critical states of thick isotropic strips and disks 
were calculated in Refs. ^3|,^4] using the appropriate 2D 
equations. The isotropic fronts (p — 0) obtained from 
Eq. (21) and depicted in Fig.Ja coincide jKith the flux 
fronts computed for such stripscia and diskscJ if one takes 
d <C a and j c = J c (0). These flux front profiles are equal 
for strips and disks since the appropriate g.tatic sheet cur- 
rents J(x) and J(r) have identical formsau described by 
Eqs. (19,20); for the disk x should be replaced by the 
radial coordinate r and a means the disk radius. Inter- 
estingly, the thin film solution can well describe features 
of thick strips and disks up to quite large aspect ratios 
d/2a « 0.2 as long as the flux front is not too close to 
the center (xq > d). 
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P = 1, q = 4 
d/ 2a = 0.04 





FIG. 6. Magnetic field lines in thin strips with aspect ratio 
d/2a = 1/25 in a perpendicular magnetic field H = 0.2j c (0)d. 
Isotropic pinning (p — 0, top) and same anisotropic pinning 
examples as in Figs. 3-5 are shown. The upper plot is to 
scale, the lower three plots are stretched along z for clar- 
ity. The depicted lines are contours of the vector potential 
with nonequidistant levels Ai oc i 2 , i = 0, 1, 2, 3, . . ., and are 
thus directed along the magnetic field. The nearly equidistant 
lines indicate a nearly constant gradient of the field magni- 
tude. The inner contour A(x,z) = coincides with the flux 
front 7. 

In a strip the sheet current has a fixed direction, and 
thus there is no rotation of the flux lines in the region 
inside 70, i.e., i\) =const in Eq. (3). On the other hand, 
since in the considered case any in-plane anisotropy of 
pinning is absent, the flux lines are not twisted in the 
region x > xq, and the solution (9-11) is valid. The 
distribution of B t = B x (x, y, z) in the strip can be found 
using Eqs. (5,11). For our model (17,18), these equations 
are solved analytically, yielding for x < xq 

B x {x,z) = 0, \z\ < z 1 (x) , 

\B x (x,z)\ = ic(0)(l +p) [\z\ - z^x)], 
z 1 {x) < \z\ < d/2, 

while in the region x > Xq we arrive at 



(22) 



\B x (x,z)\ =j c (0)\z\ 



1 + p exp 



q B z {x)d 



2\z\B c 

\z\ < d/2. 



(23) 



Note that in the anisotropic case (p ^ 0) B x is not linear 
in z in the region x > xq. Taking into account these 
equations and the fact that B z is independent of z, one 
finds the 2D distributions of B inside the thin anisotropic 
strips shown in Fig. 6. As can be seen_in the isotropic 
case our approach gives the expectedEj magnetic field 
lines except in a narrow region near the flux front, where 
this approximation breaks down (see Sec. II). The width 
of this region is less than, or of the order of d. 




FIG. 7. Magnetization curves of thin strips with 
anisotropic pinning, Eqs. (17,18). The dotted line is for 
isotropic pinning (p = 0). The loops with a hump are for 
p = 2 and the loops with a dip for p = —0.5. The dashed 
curve has half the sweep amplitude. The applied field H is 
in units of j c (0)d and M is the magnetic moment per unit 
length in units of j c (0)da 2 . Creep exponent n = 51. 



Knowing the direction of H(x,z) we can further cal- 
culate the distribution of j = j c in the region outside 
the 3D front 7 (inside 7 one has j = 0). When x < x , 
the flux lines are practically parallel to the xy plane, 
and thus j c = j c (7r/2) = j c (0)(l +p). If x > x 0l j c is 
given by the first equation (18) with t replaced by the 
ratio B z (x)/[fioj c (0)z]. Note that for p > and x > x , 
j c (x, z) is greater near the surfaces (|z| ~ d/2) than in the 
central plane of the strip [z = 0). The opposite situation 
occurs when p < 0. 



In Fig. 7 the virgin magnetization curves and magnetic 
hysteresis loops M(H) are presented for strips with vari- 
ous anisotropics of pinning. For ab plane pinning (p > 0) 
a so-called central peak appears. Thus, in high-T c super- 
conductors this often observed peakn3 in principk may 
be caused by intrinsic pinning by the CuO planes £3 Note 
that in this case the width of the peak is proportional to 
the thickness d of the sample (our field units are j c (Q)d). 
From this one may estimate the contribution of the in- 
trinsic pinning when the peak is observed. In the oppo- 
site case p < 0, i.e. when j c {9) has a minimum at — tt/2, 
a dip in \M(H)\ occurs in the region of weak fields. An 
additional dependence of j c on LBJ may then lead to a 
so-called fishtail (or peak) effecttil with a maximum of 
\M\ at some intermediate field. In Fig. 7 the zig-zag 
sweep of the applied field H(t) has constant ramp rate, 
\dH / dt\ = 1 in units where j c {0)d = a = fiQ = Eq = 1 
in the used current-voltage law E(J) = E (J/J c ) n with 
n = 51, cf. Appendix B. When the creep exponent n is 
lowered to a more realistic value n = 11 measured e.g. 
in Ref. G9, the amplitudes of the peak and dip slightly 
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decrease, and the sharp edges of M(H) at the largest and h= 
lowest H a are rounded, see also the figures in Refs. |l9|d, rj.2 




H/J = 0.2 



H/J =0.8 



FIG. 8. Stream lines of the sheet current in a thin rectan- 
gular platelet with side ratio 1.4 and isoptropic pinning. Bean 
model J c = const. Applied field H/J c = 0.2 (solid lines) and 
H/J c = 0.8 (dashed lines). Creep exponent n = 21. The 
arrows perpendicular to the stream lines and of length pro- 
portional to the sheet current, indicate the direction of the 
flux lines lying at different depths. 



V. THIN RECTANGULAR PLATES 



The critical state of rectangular platelets is qualita- 
tively similar to that of the strip. In particular, for aspect 
ratios b/a > 1.4 (see Figs. 1, 8, 9) the profiles J y (x, 0) 
and B z (x,0) along the shorter axis are practically iden- 
tical to the profiles in long strips with b/a — > oo. There 
are, however, several qualitative differences as compared 
to the magnetic behavior of long strips or circular disks: 

a) In platelets with shape different from strips or circu- 
lar disks, the stream lines of the sheet current 3(x,y) — 
{Jx, Jy) do no longer coincide with the contour lines of 
B z (x, y) and the penetrating flux fronts 70 (Fig. 1). This 
new feature is described in detail in Ref. 4 for films with 
elliptical shape. 

b) In increasing applied field H, the direction of the 
sheet current 3(x,y) changes with H at any given point 
(x, y) positioned inside the front 70 and away from the 
symmetry axes x = 0, y — 0, or edges x = ±a, y = ±b. 
In the isotropic case this rotation of J(x, y) comes to a 
halt when the flux front passes through the point. 

Both features a) and b) occur only in perpendicular 
geometry. In long rods with the same rectangular cross 
section in longitudinal field, the B z contours, J stream 
lines, and flux fronts are concentric rectangles, and inside 
the flux front both B z (x, y) and J(x, y) vanish. From the 
feature b) follows a further interesting effect: 




current stream lines 



FIG. 9. Stream lines of the sheet current J(x,y) (left) 
and contour lines of the perpendicular magnetic induction 
B z (x,y) (right) in a thin rectangular plate with out-of-plane 
anisotropic pinning and side ratio 1.4. The anisotropy j c {0) 
is described by the model Eqs. (17,18) with p — 2 and q — 4. 
The applied magnetic field is (from top to bottom) H = 0.2, 
0.4, 0.6, 1.5 in units of j c (0)d. Creep exponent n = 21. 



In thin rectangles a rotation of the penetrating flux 
lines may occur inside the front 70. In this region the flux 
lines surround a thin flat field-free core and are nearly 
parallel to the x, y plane. However, at a given point 
(x, y) the flux lines lying at different z can have different 
in-plane orientations. This rotation of flux lines is due 
to the rotation of the direction of the sheet current and 
occurs in any isotropic or anisotropic thin superconduc- 
tor with shape different from a strip or circular disk. An 
example of this rotation is shown in Fig. 8. 

The stream lines of the sheet current and the contours 
of B z (x,y) in a thin rectangular plate with strong out- 
of-plane anisotropy of pinning are depicted in Fig. 9 for 
increasing applied field H . The stream lines are contours 
of the function g{x,y) introduced in Appendix B. In all 
contour plots in Fig. 9 the same constant level spacing 
is used: Ag = 1/20 and AB Z = 1/10 in units where 
j c (0)rf = a = jiQ = 1. In our computation based on 







Ref. we used \dH/dt\ — 1 and E = 1 and n = 51 
in E = E (J/J c ) n . The shape of the depicted lo-fifps is 
very similar to their shape in the isotropic plate,tlQ i.e. 
different J C (B Z ) dependences lead to almost the same 
penetration pattern. This is so since the shape of these 
loops is determined mainly by the specimen shape: At 
small H <C J c (0), one has almost complete field expul- 
sion and the stream lines are thus independent of any 
material property, being rectangular near the edges and 
circular near the center of the rectangle. After the flux 
front has passed through a given point, J abruptly tends 
to nearly a constant, see Fig. 3 and Eq. (17). The situa- 
tion then becomes similar to the isotropic case, and the 
stream lines practically coincide with concentric rectan- 
gles with constant spacing proportional to 1/J. 

Since the current stream lines in the rectangular 
platelet do not coincide with the lines B z =const, and J 
depends on B z , the penetration pattern in the anisotropic 
case could differ-from that obtained in the isotropic 
superconductor .ErQ However, even a close look can hardly 
discover a qualitative difference in the shapes of the con- 
tours in isotropic and strongly out-of-plane anisotropic 
(Fig. 9) rectangular plates, in_stark contrast to the situa- 
tion with in-plane-anisotropy.tH A minor difference is that 
in the partly penetrated state (H = 0.4 and H = 0.6 in 
Fig. 9) the stream lines outside the flux front and near 
the axes x = and y = 0, are slightly convex towards 
the center, while for isotropic Bean pinning the stream 
lines are straight there (or are slightly concave due to 
the finite creep exponent n). As expected from the strip 
results, Fig. 3, with out-of-plane pinning the density of 
the sheet-current stream-lines near the flux front is larger 
(when p > 1) or smaller (when p < 1) as compared to 
the case of isotropic pinning {p = 0), cf. Eq. (17). The 
same statement is true for the contours of B z near the 
front. 

Figure 10 shows magnetization curves of thin rectan- 
gular plates with side ratio 1.4 and various out-of-plane 
anisotropies. These magnetization loops exhibit central 
peaks or dips caused by the anisotropy and are very sim- 
ilar to those shown in Fig. 7 for strips. Thus, one can 
conclude that the features of the loops at H ~ j c d are 
specified mainly by the anisotropy of flux-line pinning 
(i.e., by p and q) rather than by the shape of the sam- 
ple. This property should be useful for determining these 
pinning parameters from experiments. 



VI. CONCLUSIONS 

In many papers thin superconductors are considered as 
2D samples with a constant critical sheet current J c or 
with J C {B Z ) — j c (B z )d if j c depends on the magnitude 
of the magnetic induction. However, when an out-of- 
plane anisotropy of pinning exists, the current density 
also depends on the angle 8 between B (the flux-line 
direction) and the normal to the sample plane (our z 



axis), and thus the critical state problem becomes three- 
dimensional. We have shown here that the general 3D 
critical state problem for thin superconductors of arbi- 
trary shape and with arbitrary 3D anisotropy of pinning 
can be separated into a ID problem, which treats the 
current density and induction across the thickness of the 
sample, and the 2D problem of thin superconductors with 
a new induction dependence of the critical sheet current 
J C (B Z ) flowing in the plane of the superconductor. The 
new J C (B Z ) dependence is determined by the out-of-plane 
anisotropy of the critical current density j c {&) and mod- 
ifies the original B z dependence and in-plane anisotropy 
of J c , if they exist. The resulting 2D thin film problem 
with given J C (B Z ) can he. salved by standard numerical 
and analytical methods.E^To Our theory generalizes the 
approach of Ref. Uu for a circular disk. 




FIG. 10. Magnetization curves of thin rectangular plates 
with anisotropic pinning, Eqs. (17,18), side ratio a/b = 1.4, 
and n = 51. As in Fig. 7, the dotted line is for isotropic 
pinning (p — 0), the loops with a hump are for p — 2, and 
the loops with a dip for p — —0.5. The applied field H is in 
units of j c (0)d and M is the magnetic moment in units of the 
saturation moment of the rectangle, (2/3) j c (0)d& 2 (3a — b). 



Knowing the solutions of these ID and 2D problems, 
the critical state in three-dimensional thin flat supercon- 
ductors (Fig. 1) can be understood. Two qualitatively 
new features of the critical state occur in this general case 
as compared to the highly symmetric cases of isotropic 
disks or strips, which are commonly considered. First, 
when the shape of the sample differs from a strip or a 
disk then a rotation of flux lines (differently oriented 
in different planes z = const) may occur in the region 
where B z = 0, i.e. where the flux lines surround a non- 
penetrated flat core and are nearly parallel to the sample 
plane since the sample is thin. Second, if both out-of- 
plane and in-plane anisotropies of pinning exist, the flux 
lines even in thin superconductors are not only curved 
but also twisted in the region where B z ^ 0, i.e. outside 
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the front line 70 in Fig. 1. Twist means that the in-plane 
component of the direction of the flux lines at a given 
point x, y changes with the depth z. 

Using the above approach, the critical states of 
strips and rectangular platelets with only out-of-plane 
anisotropy are analyzed. We consider a rather general 
analytical model (18) for this anisotropy, which can sim- 
ulate both the intrinsic pinning by CuO planes and pin- 
ning by extended defects aligned with the c axis in high- 
T c superconductors. We obtain the shape of the flux 
front and the two-dimensional distributions of the mag- 
netic induction and current density for the anisotropic 
strip. In the case of intrinsic pinning, the flux front has 
a wide flat part near the central plane, Fig. 5; the cur- 
rent density peaks near this flat front, and the induction 
profile is very .steep there, see Fig. 3. In fact, a so-called 
current stringcZrEj occurs in this case, i.e., an additional 
current flows along this steep front as compared to the 
isotropic situation. In Refs. g7 28 the string was consid- 
ered in connection with a finite lower critical field H c \. 
However, in the considered case the string is not due to 
a finite iJ cl but results from the anisotropy of pinning, 
and the jump of B z at the flux front considerably exceeds 
H c i since H cr = j c d/2 3> H c \ was assumed above. This 
type of out-of-plane anisotropy leads to a central peak in 
the magnetization curve, Fig. 7. 

A different situation occurs when the 9 dependence of 
j c has a minimum at 9 = tt/2. In this case the flux front 
has the shape of a wedge with a rounded point that be- 
comes sharper with increasing anisotropy. As can be seen 
in Figs. 3b, c the current density near the front line 70 now 
has no peak and no sharp step (as for isotropic pinning) 
but it decreases monotonically and has only a rounded 
step and an inflection point with vertical slope at this 
front line. The induction B z now vanishes less steeply 
than in the isotropic case, and in the limit of large such 
anisotropy the induction profile near the front decreases 
almost linearly. The magnetization loop for this type of 
anisotropy has a central dip, see Fig. 7. 

All these features remain qualitatively the same when 
a smaller creep exponent n is chosen in E oc J", Eq. (B6), 
or when the superconductor has a rectangular shape. In 
particular, with smaller n the cusp in J(x) remains sharp 
and the profile of B z (x) remains steep, see Fig. 3d with 
n = 11. As shown in Figs. 7 and 10, the magnetiza- 
tion loops of thin rectangular plates with out-of-plane 
anisotropy exhibit the same central peak or dip as thin 
strips. 

For any given dependence j c (B), the equations of this 
paper enable one to compute various characteristics of 
the critical state in thin superconducting samples of re- 
alistic shapes. With this theory at hand, experimental in- 
vestigation of flux-density profiles, of the H dependence 
of the penetrating flux front (e.g., by magneto-optics), 
and of magnetization loops, will yield information not 
only on the strengh but also anisotropy of flux-line pin- 
ning in superconductors. 
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APPENDIX A: ELECTRIC FIELD IN 
ANISOTROPIC SUPERCONDUCTORS 

In a uniform anisotropic superconductor the direction 
of the electric field E generated by vortex motion in gen- 
eral will not coincide with the direction of the current. 
To study this effect, we start with the situation when a 
current with the density j flows in the plane normal to 
the magnetic induction B (below we call this plane the N 
plane). The creep activation barrier U{B , tp, j / jo(B , tp)] 
is assumed to vanish at j = jo(B,tp) [this is the defini- 
tion of the factor jo(B, tp)] . Here the angle tp specifies the 
direction of the current in the N plane, j cx (cos tp, sin tp), 
and U[B, <p,j/jo(B, tp)] is the value of the barrier which 
a flux bundle has to overcome to hop in the direction of 
the Lorentz force, i.e., along the vector [j x B]. At given 
temperature T, the expression exp{— U[B, (, j cos(C — 
tp)/jo(B, C)]/T} gives the probability for a flux-bundle 
to jump in the direction normal to any given in-plane 
vector (cos C, sin C). Hence, the electric field E is ob- 
tained by averaging this expression over all angles C in 
the interval tp — tt/2 < C < p + tt/2 (we disregard flux- 
bundle jumps against the Lorentz force). If the effective 
depth of the flux-pinning well [dU(B,tp,x)/dx] x =i con- 
siderably exceeds T, the direction of E in the N plane, 
(cos Co 7 sin Co), can be found as the position of the mini- 
mum of U[B, Q,j cos(( — tp)/jo(B, ()] with respect to vari- 
ations of C- Since, according to the definition of jo(B,Q, 
one has dU(B,£,l)/d( = 0, one obtains the following 
formula for Co in the critical state: 



tan(Co - tp) 



fo(BXo) 

io(£,Co) ' 



where j' (B, C) = djo(B, C)/^C- As for the critical current 
density j c (B, tp), we have 

Jc (B,p) = io( .f' Co) . = [j„(B,Co) 2 +J„(3,Co) 2 ] 1/2 • 
cos(Co-^) 

It is also useful to express the angle Co m terms of 
j c (B,p), i.e., the quantity determined in experiments, 

S = ( -p = - arctan ( J '^f ' ^ \ , (Al) 



jc(B,tp)J ' 



where j' c (B,tp) = dj c (B,p)/dp. The above formulas 
have been derived under the assumption that cos(C — 
tp)/jo(B z , C) has only one maximum in the interval tp — 
tt/2 < C < p + n/2. This situation does not always occur 
in twinned crystals. In such crystals the angle Co m &y 
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have a fixed value in some interval of angles ip, and thus 
j c (<p) — Jc(Co)/ cos(Co — (p) inside the interval. However, 
note that Eq. (Al) is still valid in this case. 

We now present the formulas for the case considered 
in Sec. II: critical currents flow in the xy plane with 
jc = {jc cos ip, j c shx(p), and a flux-line element has the 
direction (sin#cos'!/>, sin 9 sin ip 1 cos (9). Since in the ab- 
sence of flux cutting only the component of j c normal to 
B is essential, the critical current density in the N plane, 
jj" , has the form 

£(B„Bt,<p-il>,il>) = n x (cp)j c (B z ,B t ,<p-iP,^) , (A2) 

where n±((p) = [1 — cos 2 (tp~ ib) sin 2 O] 1 ^ 2 . It follows from 
Eq. (Al) that the angle 6 between E and j,f is given by 
the formula 



with 



tan 5 = — 



d[lnj^(B z ,B t ^-^^)} n 2 ± (<p) 



(A3) 



dip cos 9 

The unit vector e directed along E has the components: 
1 



1 



[cos £ — sin 9 cos ip cos(£ — ip)} 



-v — if\ I sm £ _ sin2 sin ^ cos (£ _ VO] > (A4) 



MO' 



1 



• sin 9 cos 6* cos(£ — ip) 



n±(0 

where the angle £ is defined by the expression 

, . „tan(o9 — ip) + cos tan 5 

tan £ - V) = cos 7 . (A5 

cos 9 — tan a tan(</? — V) 

In other words, e is the normalized projection of the vec- 
tor (cos£, sin£) onto the N plane. This explains the 
meaning of the angle £. 

It is important to keep in mind one special case. If an 
in-plane anisotropy of flux-line pinning is absent, it fol- 
lows from symmetry considerations that at fixed B z , B t , 
ip the derivatives dj^r/dip and dj c /dip are equal to zero 
when ip — ip — ±7r/2. Thus, in this case one has 5 = 0, 
£ = ip, and the directions of E, j c and ]jr coincide. 



APPENDIX B: DYNAMIC EQUATIONS 

It is convenient to express the sheet current in the film 
through a scalar function g(x, y) = g(r) as 

J(r) = V x zg(r) = — z x V.g(r) . 

This substitution guarantees that V • J = and that 
the current flows along the specimen boundary T if one 
puts g(r) = const = there (the lines g(r) = const co- 
incide with the current-stream lines). Then, Eq. (14) is 
transformed as follows ,E3 

^^ + C ( r). 9 (r)-/^-|W rf V, (Bl) 
Mo Js 



2tt 



R(cp) 



(B2) 



The integral in Eq. (Bl) is taken in the Cauchy sense. 
In formulas (B2) the integration area S is the entire a;, y 
plane with exclusion of the sample area S; R{4>) = |s — r|; 
s defines the position of a point on the boundary T; and 
(p is the angle of the vector s — r relative to an arbitrary 
fixed direction. Eq. (Bl) admits inversion, 



g(r)= f Kiv^^B^-H^r' . 
Js 



(B3) 



where the kernel K(r, ri) is found from the following 
equation, 



C(r) 

K(t',t 1 )-K(t,t 1 ) j 2 , 

47Ti?3 



} 



(B4) 



which can be solved by iteration starting with K (r, r-J = 
0. 

We now differentiate Eq. (B3) with respect to time 
t, take into account the Maxwell equation, dB z /dt = 
i • [V x E], and the current-voltage law,E2l 

E = E(J,£ Z ), 

which expresses the electric field E through the sheet 
current J. Then, we arrive at an equation for <?(r), 



dg(r,t) 
dt 



I K(t, r') 
Js 



is 

-[zxV].E[-zxVj(r',(),5 z (r')] 



dH(t)\ M l 



dt J 



d\. (B5) 



Equations (Bl) and (B5) allow one to find B z (r, t), g(r, t) 
[and thus J(r, t)] if the prehistory H(t) and the current- 
voltage law are given. 

The static critical state can be calculated by a dynamic 
approach using any model law E which has a sharp bend 
at J = J c . The specific form of this dependence is irrel- 
evant; e.g. we may use the power law 



(B6) 



E — En ( — 



and take the limit n 3> 1. That the static critical state 
is indeed reached in the limit n — > oo was shown rigor- 
ously in Ref. 31 The function J(B Zl ipo) resulting from 
Eq. (12) or from the appropriate generalization of this 
equation should be used as the critical sheet current J c 
in Eq. (B6) for currents with direction (cos^OjSin^o)- 
The direction of the constant vector E in general does 
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not coincide with J if there is anisotropy of flux-line pin- 
ning in the xy plane (see Appendix A). At given tpo, the 
angle £ defining the direction of E oc (cos £o, sin £o) is 
determined by the formula 



tan(£ - <A)) 



d[]nJ(B z ,<p )] 



(B7) 



Thus, to describe the critical state of an infinitely thin 
superconductor, it is sufficient to solve the dynamical 
Eqs. (B1,B5,B6) with the initial condition g(x,y)\t=o — 
0. This dynamical approach appears to be more conve- 
nient than solving the static Eqs. (13,14) with an initially 
unknown flux- front 70. 
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